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The impact of local symmetries on the transport properties of one-dimensional aperiodic arrays of 
potential barriers is explored. We show that weak commutation of combined local parity operations 
with the Hamiltonian induces states which either possess definite local parity and zero current or are 
totally transmitting with their probability density following the local symmetries of the potential. 
Our findings are illustrated through a construction principle for barrier arrays with prescribed 
transport features. 

PACS numbers: 03.65.-w, 73.40.Gk, 03.65.Xp, 73.63.Hs 



Introduction. — Symmetries set the grounds on which 
physical systems are treated, their ubiquitous presence 
on all scales dictating the form of the developed theory 
as well as experimental observation and analysis. Global 
spatial symmetry, though, in most cases pertains exactly 
to structurally simple systems and idealized models, while 
its absence in realistic situations often implies their sta- 
tistical treatment To quantify situations of approxi- 
mate symmetry, measures of symmetry [2, 3] have been 
proposed, which reflect the degree to which it is fulfilled 
under specific operations. In the generic case, symme- 
try of a system under spatial transformations is globally 
broken, but retained on a local scale [4], such that the 
associated invariance of physical properties affects its be- 
havior beyond a mean description [5]. Local symmetries 
can be present by design, as in mesoscopic transport de- 
vices [6] or other artificial experimental setups [7], but 
also naturally inherent in structurally complex systems, 
like large molecules 0, quasi-crystals [9] or even disor- 
dered matter [lo|. In view of the intensive exploration 
of such systems, the need for a rigorous theoretical treat- 
ment which addresses local regularities, but incorporates 
the global composite structure, becomes evident. 

In this letter we make a step in this direction, inves- 
tigating the impact of local symmetries in the transport 
properties of globally non-symmetric systems. Consider- 
ing scattering in one dimension (ID), any spatial sym- 
metry transformation can be reduced to combinations of 
coordinate inversions within subintervals of configuration 
space, which we refer to as local parity (LP) transforma- 
tions. 

Utilizing non-uniform arrays of arbitrary potential bar- 
riers, we show how LP is related to the transport prop- 
erties of stationary scattering states with symmetric or 
asymmetric asymptotic conditions (SAC or AAC): un- 
der conservation of LP, the total or reflected probabil- 
ity current vanishes, respectively. The latter corresponds 
to totally transmitting resonances (TTRs) in the non- 
symmetric system which we classify with respect to the 
symmetries of the scattering state profile. Applying the 
concept of LP explicitly to the analytically tractable case 
of rectangular barriers, we demonstrate a construction 
principle for aperiodic devices with prescribed transport 
properties. 



Local parity. — Let us first introduce the LP operation, 
which performs the usual parity transform in a finite sub- 
domain V of configuration space (of the x-axis in ID) 
and acts in the remaining part, up to a sign, as the iden- 
tity operator, in order to maintain spectral equivalence 
to global parity. The single LP operator 77^ = IJ^ ± is 
parametrized by the location of its inversion point a and 
the size L of P, in terms of which it is defined by its 
action on an arbitrary function f{x) as 

f(2a-x) 0{L/2- \x-a\) 

±f{x) 0{\x-a\-L/2), (1) 

with eigenvalues X± = ±1, where O denotes the Heav- 
iside step function. Two different LP transforms il^S 
77^^ commute if the associated domains do not overlap, 
Pi n P2 = 0- Subsequent application of N = A/'+ + N- 
non-overlapping single LP transforms, where N± is the 
number of transforms with parity = ±1, thus corre- 
sponds to a total LP operator = H^^-i with 



depending on the number of odd 



eigenvalue A = (—1)^ 
transformations. 

The link between local symmetry and transport prop- 
erties is provided by the commutation of the Hamilto- 
nian H = P + F (using units h = 2m = 1) with 
71^°*, and thereby with each TI^J". For a potential 
V{x) which is mirror symmetric within the subdomains 
= [o^n — Ln/2^an-\-Ln/2]^ the actiou of the associated 
commutators reads 



[H,n^:M{x) 



{f{2an-x)Tf{x))A\x) 
-4 {f\2an-x)±f\x))A{x) (2) 

where A(x) = S{x — an — Ln/2) — S{x — an-\-Ln/2) and the 
prime denotes ^-differentiation. Note that, although V 
clearly commutes with fl^'^ , the kinetic term in H leads 
to non- vanishing boundary terms, ordered in degree of 
singularity. The commutation is then manifest in a weak 
sense, that is, subsequent action of H and each 77^ on 
f{x) is independent of their order, only if the boundary 
conditions 



f{an-Ln/2) = ±f{an^Ln/2) ; 
f{an-Ln/2) = Tf{an^Lj2) , 
are fulfilled for each n = 1, 2, A^. 



(3) 
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Scattering. — We now explore the consequences of local 
symmetries in ID scattering, based on the derived condi- 
tional commutation of H with 77^°* . The generic setup is 
defined by the real static potential 



V{x) 



Ns 



e{L„/2 



(4) 



describing an array of Ns non-overlapping scatterers of 
widths centered at and symmetric about 
{Vn{2an — x) = Vn{x)). SAC or AAC are imposed on 
the stationary scattering state {ipk) with plane waves 
^±zkx _ (^x\^) of energy E = incident on the left and 
right, or only on the left, of the scattering region, respec- 
tively {k > 0). The corresponding asymptotic ingoing 
amplitudes (}) or (J) of \ijjk) in |±)-space are scattered 
off the potential into the outgoing ones by the S'-matrix 



r t 
t r 



i'd 



-y/Re- 



i5 



(5) 



which is unitary by momentum (probability current) con- 
servation and symmetric due to the time-reversal invari- 
ance of our system. T = |tp = 1 — |rp = 1 — R is the 
transmission probability, which is measurable for a setup 
under AAC. 

Zero- current states. — Consider first the scattering 
setup under SAC, with wave function i^^{x) = e*^^ -\-{t-\- 
r)e~*^^ on the left {x < ai — Li/2) and {x) = e~*^^ + 
{t + f)e*^^ on the right {x > ^ of the array. 

A realization of this could be thought of as a split parti- 
cle beam brought to interference at the potential region. 
We now imply the conservation of LP in the subdomains 
= — Qfn + ^] (n = 1, 2, Ns) of the scatter- 



ers and the subdomains 



\OLn 



2 -> ■ 



(n = 
-1 for 



1, 2, Ns — ^) between them, with A^, = +1 or 
even or odd transformations. This imposes the conditions 
dU with f{x) = ipkix)^ N = 2Ns — 1, on the total wave 
function at adjacent subdomain boundaries, and thereby 
at the left and right boundary of the potential by induc- 
tion: 

Mo^i - Li/2) = {-l)^-Mo^Ns + Lns/2) ; 
V^Uc^i - Li/2) = -(-l)^-V^Uc^iv, + Lns/2) , (6) 

tjjk here being a common eigenfunction of H and IJ^^. 

These in turn lead to the following conditions for the 
energies of the LP eigenfunctions and for the form of the 
S'-matrix: 



=(-1) 



(7) 



^r-\-(^Ns~^ — ^ )/2is the center of the ar- 



where Xc = (o^i 2 ' ^^^s ' 2 
ray. This means that, at momenta ki^^ = n7r/2xc, n G N, 
for a given location Xc, the locally symmetric array be- 
haves as a globally symmetric one with respect to its 
asymptotic transport properties. In particular for the 
setup under SAC, at these energies the probability den- 
sity current jk = 2k{l — \t-\-r\) = —2k{l — \t-\-r\) vanishes 
due to the unitarity condition S = I with r = f (7^ in 
general). This can also be understood from the fact that 



10- ;Re(V^fcJ/fcc Im(^fcJ//ce 



PkjEr 



zero-current state 
(SAC) 
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FIG. 1. Transmission spectrum of the simplest LP-symmetric 
barrier setup, with scaled ipkd^) and pkr^{x) shown for a zero- 
current state {Ec = 8.70) under SAC and a TTR {Er = 4.60) 
under AAC, respectively. 



r = f yields a phase difference 5 = 7r/2 between t and r, 
as seen from Eq. (j5j), so that |t + r| = vl^P^^kF = 1- 
The zero-current states for SAC connect the concept of 
LP to a transport observable, in a situation where the 
transmission coefficient T is ambiguous due to the indis- 
tinguishability of reflected and transmitted wave ampli- 
tudes in the asymptotic regions. LP conservation, and 
thus j/e, is sensitive to the relative phases of the ingo- 
ing amplitudes (or, equivalent ly, to the positioning of the 
array) at a flxed energy, thus implementing a ID inter- 
ferometric setup based on probability current. Examples 
of LP-conserving zero-current states are shown in Figs. [1] 
and [3] for the case of rectangular barrier arrays (see be- 
low). 

Totally transmitting resonances. — Keeping the LP- 
sensitive current jk as an observable, we now turn to the 
case of AAC, which are most common in transmission 
measurements. With incoming amplitudes (J) in 
space, that is, '^^(x) = e*^^ + j,^-%kx '^^(x) = te*^^, 
the symmetry in the scattering problem is broken ex- 
plicitly, so that tl^k{x) will clearly not conserve LP un- 
der AAC. On the other hand, the (spatially independent) 
current is now directly connected to the probability den- 
sity of the transmitted wave: jk = ^IV^^P = kT. We 
therefore seek a remnant of the conserved LP in the spa- 
tially resolved probability density pk{x) = |V^/c(x)p of the 
scattered particle with AAC. To this aim, we repeat the 
arguments leading to conditions ([3]), setting f{x) = pk{x) 
and restricting the LP operation to even transformations 
(An = +1 or N- = 0), since pk{x) > by deflnition, 
which leads to the corresponding conditions by induction 
through the array: 



Pki^i - Li/2) = pk{c^Ns + Lns/2) ; 



(8) 



To see the consequences of the even LP of pk{x) within 
each of the array subdomains, we impose the above con- 
ditions on the (squared) amplitudes connected by the S- 
matrix. This yields 



(a) 







(b) 



„2ife(ai-Li/2) 



(9) 



corresponding to total transmission (T = 1) or total re- 
flection {R = 1), respectively. Concentrating on the for- 
mer case pjj], we thus conclude the following: The scat- 
tering state under AAC exhibits a TTR if and only if its 
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probability density is non-zero everywhere and LP sym- 
metric with respect to the local symmetry axes of the po- 
tential. We will demonstrate examples of TTRs in rect- 
angular barrier arrays without global parity symmetry 
below. 

This result generalizes the relation between resonant 
transmission and symmetric probability density, from 
globally to locally symmetric potentials, in a conceptu- 
ally transparent way. For a single symmetric scatterer, 
an isolated [HI resonance at momentum k = kr can 
be shown to have symmetric pk{x) and, therefore, be 
totally transmitting, by expanding the propagator of 
the wave function into resonant global parity eigenstates 
[l3| . This confirms the intuitive symmetry argument, 
that the (observable) probability density of a totally 
transmitting state cannot reveal the direction of inci- 
dence in a reflection symmetric potential. The symmetry 
of resonant states is not directly evident, though, if 
different such (symmetric) scatterers are connected into 
a non-symmetric array. Only when treated separately, 
each symmetric scatterer could be argued to transmit 
resonantly at the same kr with globally symmetric pk{x). 
This would then remain unaffected by implementation 
into an array, since only phase-shifted plane waves 
propagate in the intervening potential- free regions. 
Nevertheless, this argument does not carry over to the 
case of SAC, where counterpropagating waves interfere 
between the scatterers. Moreover, at a TTR not all 
local symmetries of the potential are necessarily obeyed 
by the resonant probability density, as we will see in 
the following. The proposed concept of LP treats the 
non-symmetric system globally while addressing its 
local symmetries in arbitrary subdomains, and applies 
uniformly to scattering under SAC or A AC. 

Rectangular harrier arrays. — With the general consid- 
erations above applying to arbitrary (smooth) potentials, 
we now demonstrate the LP induced zero-current and 
TTR states in the analytically solvable case of rectangu- 
lar barrier arrays, that is, Vn{x) = V^, n = 1,2,...,A^5, 
in Eq. (|4]). The potential strength V^, width Ln and lo- 
cation an of the n-th barrier are freely adjustable. The 
barriers themselves and the regions between them define 
the smallest possible subdomains with local symmetry 
and mutually different potential strengths. Relying on 
regularities at the level of barriers, though, local symme- 
tries can also be realized on larger scales, which is gener- 
ally the prevalent case in natural complex structures or 
fabricated devices. 

We therefore group the Ns barriers into < Ns 
reflection symmetric resonators IZi (I = 1,2,..., N^), sep- 
arated by regions {l = 1,2,..., A^t^ — 1) of zero poten- 
tial. The resonators are of three types with respect to 
the TTRs they support: (I) a single barrier, supporting 
above-barrier resonances (ABRs), (II) a homogeneous ar- 
ray of identical equidistant barriers, supporting tunneling 
resonances and ABRs, or (III) an inhomogeneous barrier 
array, supporting isolated TTRs for appropriate combi- 
nations of barrier strengths and widths |12|] . 

For any resonator decomposition, the |±)-amplitudes 
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FIG. 2. Setup of two large resonators IZi (with m = 5 barri- 
ers) and 7^2 (^2 = 9), exhibiting two close (tunneling) TTRs 
at El = 7.00 and Ei = 7.46, with corresponding probability 
densities shown for A AC. 

along the array are determined by matching the 
wave function at the 2Ns interfaces, and connected 
by the single-barrier unimodular transfer matrices 
Mn{k; V^^Ln)' The total transfer matrix is then advan- 
tageously given by the (ordered) matrix product 

^ = r* J* = n ^-(^5 V^, Ln) (10) 

^ ^ n=l 

connecting the amplitudes of ?/^^ on the right of the array 
to those of on the left, where the matrix elements are 
related to the reflection and transmission amplitudes of 
the ^'-matrix hy w = (t*)"^ and z = -r*(r)"^. 

Conditions r = r and r = for zero current k = kc 
(SAC) and TTRs at /c = /c^ (AAC) are now expressed as 
an imaginary or zero transfer matrix element respec- 
tively: Re{z) = (SAC) or z = (AAC). In practice, 
though, conditions (|3]), with f{x) = ^l^kix) or pk{x), are 
imposed directly at the boundaries of the subdomains Vi 
with LP, which contain the resonators TZi, and of the Vf 
between them. TTRs then resonate within and between 
the resonators for a specific kr. The decomposition of the 
array into resonators in general changes for different res- 
onant energies, according to the local symmetries of the 
corresponding probability density p/c^(x). If the system 
resonates within smaller constituents of a resonator, then 
we refer to the resonance as reducible. The invariance of a 
given resonant energy under interchange or translation of 
such constituents implies their independence within the 
total system [14], in terms of transmission. In this sense 
also resonators are independent constituents of the array 
at TTRs. 

The fulfillment of LP through the conditions imposed 
on the transfer matrix suggests a construction principle 
for the design of aperiodic arrays with TTRs at prescribed 
energies, the basic steps being as follows: First, we calcu- 
late the element zi{k; {V^^^^^}, {^f^^}) of ^^e transfer ma- 
trix Ml for each resonator TZi containing ni barriers of 
strengths V^''^^ and widths Lf\ i = l,2,...,n/. We then 
impose the conditions (|3|) for even LP of f{x) = pk{x) 
within each resonator, leading to a set of 2A^7^ real alge- 
braic equations for given resonant momentum kr, 

zi{kr;{V°^'^},{Lh) = 0, I = 1,2,..., Nn. (11) 




FIG. 3. Multi-barrier array decomposed into Nn = 5 res- 
onators (solid lines), exhibiting a zero-current state (with 
shown Re(V^/eJ for SAC) at Ec = kl = 7.354. If 7^l is re- 
moved (dashed lines), the state becomes a TTR (with shown 
pkc for AAC) at the same energy. 
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FIG. 4. Setup of Ns = 12 rectangular barriers, two of which 
are deformed in width and strength. The array is decomposed 



into N. 



(1) 
7^ 



5 and N. 



(2) 

7^ 



12 resonators with respect to sup- 



ported irreducible TTRs at energies Ei — 5.00 and E2 — 8.47. 



The 2 Yld^ setup parameters can be uniquely deter- 
mined only for = 1, that is, each resonator being of 
type I. However, since generally ni > 1, it is in prin- 
ciple possible to obtain TTRs at multiple energies, or, 
conversely, construct different aperiodic arrays which res- 
onate at the same kr. Note that, for the case of SAC, also 
the intervening potential- free regions need to be adjusted 
to conserve LP. 

Representative setups. — Let us now discuss some ex- 
amples which illustrate the above concepts and their im- 
plementation more explicitly, in order of increasing res- 
onator number A^t^. Fig. [T] demonstrates the occurrence 
of a zero-current state and an irreducible TTR in the 
simplest case which breaks global parity, an asymmetric 
double barrier setup. As we see, the zero-current state 
has even LP in the barriers and odd LP between them. 
T{E) displays a multi-ABR structure, but only the de- 
picted peak corresponds to a TTR state (see inset). 

Fig. [2] shows again a setup of N-jz = 2 resonators, but 
of a more complex structure, supporting two energetically 
close TTRs. (at E2) is irreducible in both resonators, 
while pk^ (at £^1) is irreducible in IZi but reducible in 
7^2, where it varies only within the barrier regions (with 
flat p{x) indicating a forward propagating wave only). We 
point out that the setup remains resonant for an arbitrary 
combination of IZi and 7^2, with repeated corresponding 
patterns in pk{x)^ which thus resembles an extended sys- 
tem with structural complexity. 

In Fig. [3] a zero-current state becomes a TTR at the 
same energy by removing the first resonator IZi from a 
N-jz = 5 array. The TTR is reducible within IZ3 and 
irreducible in 7^2, 7^4, 7^5, while more localized within IZ4. 
Notice that this TTR is independent of the width of the 
central barrier in 7^5, within which the wave propagates 
only forwardly with constant pknix) > 1. This might be 
utilized for the flexible design of efficiently transmitting 
non-symmetric devices. 

In Fig. |4] we investigate a (finite) periodic array with 
two defects in the form of alternate lattice cells. Without 
the defects, the transmission properties of the uniform ar- 
ray are determined by its unit cell , which also defines 
the scale of local symmetry. Due to coupling of the de- 
generate resonant levels of adjacent identical resonators. 



a uniform A/'-scatterer array exhibits {N — l)-fold split 
TTRs [15|, ll6|, which saturate into transmission bands 
for larger N. As we see, the presence of aperiodicity [l^ 
distorts the precursors of the energy bands [l|, [lij and 
lowers the resonances in T{E) from unity because of the 
induced asymmetry [6|. Nevertheless, the decomposition 
into resonators (particularly of type III) containing mul- 
tiple barriers reveals the possibility of TTRs, as explained 
above, owing to the locally symmetric pk^ of (irreducible) 
resonant states. This implies that identification of local 
symmetries on scales larger than the lattice cell provide 
a key for the description of structurally complex systems. 
Moreover, different LP axes (i.e. resonator decomposi- 
tions) for the same setup correspond to different TTR 
levels; local symmetry considerations thus prove to be of 
fundamental importance in accessing and understanding 
the properties of no n- uniform systems. 

Conclusions. — In conclusion, we have introduced the 
concept of local parity and revealed its impact on the 
transport properties of aperiodic ID arrays of reflection 
symmetric scatterers. LP was shown to underly the 
emergence of zero-current states and totally transmitting 
resonances under imposed symmetric and asymmetric 
asymptotic conditions, respectively. Based on these 
results, a construction principle was developed for the 
flexible realization of mesoscopic devices with designed 
transport properties relying on their LP symmetries. 
The decomposition of globally non-symmetric arrays 
into different locally symmetric resonator units relates 
total transmission to symmetries of the probability 
density. This in turn demonstrates the importance of 
considering local order on different scales to understand 
the behavior of structurally complex systems in general. 
On the other hand, the reducibility of resonant states 
into subdomains of defined resonators links to the 
concept of independence among constituents of extended 
systems on symmetry grounds. The generalization of 
our approach to higher dimensions and different kinds 
of local symmetry transformations could provide a novel 
context for the analysis of complex systems based on 
fundamental principles. 
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